Universal properties of hard-core bosons confined on one-dimensional lattices 
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Based on an exact treatment of hard-core bosons confined on one-dimensional lattices, we obtain 
the large distance behavior of the one-particle density matrix, and show how it determines the 
occupation of the lowest natural orbital in the thermodynamic limit. We also study the occupation 
Xr) of the natural orbitals for large-r; at low densities. Both quantities show universal behavior 
independently of the confining potential. Finite-size corrections and the momentum distribution 
function for finite systems are also analyzed. 
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Trapped atomic gases at very low temperatures be- 
came in the past years a center of attention in atomic and 
condensed matter physics. Particularly interesting is the 
case where the dynamics of the system is restricted to 
one-dimension (ID) due to a strong transversal confine- 
ment. It has been shown recently [ij that in regimes of 
large positive scattering length, low densities, and low 
temperatures, a quasi-lD gas of bosons behaves as a 
gas of impenetrable particles, i.e., as hard-core bosons 
(HCB). Ultracold Bose gases in ID have been realized 
experimentally |3l , so that it is expected that soon it will 
be possible to make the ID HCB gas a physical reality. 

The ID gas of HCB was first introduced theoretically 
by Girardeau j3| , who also established its exact mapping 
to a gas of noninteracting spinless fermions. Since then, it 
remained a subject of recurring attention, and a number 
of exact results were obtained for the momentum distri- 
bution function (MDF) n(fc) and the one-particle density 
matrix (OPDM) p{x) in the homogeneous 0, 0| and the 
periodic case. It was shown that n{k) ~ |fc|~^/^ for 
fc 0, and that such a singularity arises due to the 
asymptotic behavior p[x) ~ Icc]"-^/^ for large cc 0, 0. 
Similar results can be obtained using the hydrodynamic 
approximation (bosonization) 0. 

Recently, the attention has been concentrated on the 
ground state properties of the ID gas of HCB confined by 
a harmonic potential 0, as a model for the experiments. 
It was found that the occupation of the lowest natural 
orbital (NO) of the trapped system scales as \/Nb, where 
Nb is the total number of particles, as in the homoge- 
neous case (the NO are defined as the eigenstates of the 
OPDM 0). The results in ID are in contrast with higher 
dimensional systems, where a Bose-Einstein condensate 
(BEC) was proved to exist ^-^d complement that 

proof by showing that in ID only a quasicondensate is 
possible. The introduction of an optical lattice opens 
new possibilities to engineer strongly interacting states 
in the trapped gases [ij. Unfortunately, much less is 
known in the lattice case with confinement. 

We present here an exact study of trapped HCB on a 
lattice. By means of the Jordan- Wigner transformation 
[T^ we calculate exactly the Green's function for large 



systems (up to 10"* lattice sites). We find that the OPDM 
Pij ~ x^^/^ [x ~\ Xi — Xj I) for large x, irrespective of 
the confining potential chosen, even when portions of the 
system reach occupation = 1, such that coherence is 
lost there. The power law above is shown to determine 
the scaling of the occupation of the lowest NO in the 
thermodynamic limit (TL). In addition we find a power- 
law decay of the NO occupations (A,,) for large- 77 at low 
densities, and show that its exponent is also universal. 

We consider HCB on a lattice with a confining poten- 
tial with the following Hamiltonian 

Hhcb = -i ^(felfo.+i + H.c.) + 14 ^ xf n„ (1) 

i i 

with the addition of the on-site constraints b]'^ = = 0, 
l^j, 5j| = 1. The creation and annihilation operators for 

the HCB are given by bj and bi respectively, rii = blb^ 
is the particle number operator, t is the hopping param- 
eter and the last term in Eq. ^ describes an arbitrary 
confining potential, with power a and strength Va. 
The Jordan- Wigner transformation. 



(2) 
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maps the HCB Hamiltonian into the one of noninteract- 
ing fermions 

Hf = -t + H.c.) + K ^ x>{, (3) 



with fl and being creation and annihilation operators 
for spinless fermions, and n{ = fjfi- 

The Green's function for the HCB can be expressed 
using Eq. Q as 

G^, = i'f'icBhbpHCB) 

= (*^ine"^^^'^/jin^""^'^^i*^)- (4) 

/3=1 7=1 
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\^hcb) is the ground state for the HCB and is 
the ground state for the noninteracting fermions. The 
latter is a Slater determinant, i.e., a product of single 
particle states 1*^) ^ jj^^^ ^^^^ p^^/t |o), with Nf 
the number of fermions {Nf = Ni,), N the number of 
lattice sites and P is the matrix of the components of 



It is easy to see that the action of 117=1 ^ 



on the fermionic ground state in Eq. generates only a 
change of sign on the elements Pa^s for f < J ~ Ij and the 
further creation of a particle at site j implies the addition 



of one column to P with the element P, 



jNf + l 



1 and all 



the others equal to zero [the same applies to the action 

of n^=\ e"^f'^"fr on the left of Eq. Then the HCB 
Green's function can be calculated exactly as 

Nf + l N Nf + 1 N 

= (oi n E p'pifp n E ^7^4 io) 



(T — 1 7 — 1 



det 



i'A 



i'B 



(5) 



A B 

where P' and P' arc the new matrices obtained from 
P when the required signs are changed and the new 
columns added. A p roof of the last step in Eq. (jSJ can 
be found in Ref. [ijl. The evaluation of Ga is done nu- 



merically and the OPDM is given by pij = {bib 
Gij + (1 - 2Gu)- The NO (0^) can be determined 
by the eigenvalue equation J2jLi Pij 'I'] — with 
being their occupations. 

We focus next on the large- a; behavior of the OPDM. 
For the periodic case {Va = 0) we obtain that for any 
density p = Nt/N 7^ 0, 1 the OPDM decays as a power 
law pij ~ Ap/y^x/a for large x (Fig. QJ, where Ap de- 
pends only on the density (a is the lattice constant). This 
behavior was found before by means of exact analyti- 
cal treatments j^. In the presence of a confining poten- 
tial, the case relevant for the experiments with ultracold 
atoms, the situation is more complicated since the system 
loses translational invariance and no analytical results 
are available. We first analyze the case where < 1 all 
over the system. We find, remarkably, that in this case 
the OPDM decays as a power law pij ^ A'^\x/a\^^^^ for 
large x, i.e., independently of the local changes of the 



density. (They become relevant only when ni 
A 



0.) 



p depends on the characteristic density of the system 
p = Niitt/C and the power a of the confining potential. 
C = iVa/t)~^^°' is a length scale of the trap in the pres- 
ence of the lattice [3. Moreover, the exponent of the 
OPDM power-law decay does not depend on the power 
a of the confining potential, i.e., it is universal (Fig.^. 

A drastic difference between the continuous case and 
the one with a lattice is the possibility to build up regions 
with densities = 1, so that such sites are not any more 
coupled coherently to the rest. Once such regions appear, 
many NO become occupied with A,, = 1, and all the other 
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FIG. 1: OPDM vs x/a {x =| Xi — Xj |) for a periodic system 
with p — 9.1 X 10~^, Nt—91 (dashed-dotted line), harmonic 
traps {a = 2) with p = 4.5 x 10"^, iV6 = 100 (dashed line) and 
p = 2.7, A^(,=501 (thick continuous line, a = 1 region is 
present), and a trap with power a = 8, p = 7.6 X 10"*, iVi,=ll 
(dotted line). In the trapped cases the abrupt reduction of pij 
occurs when Uj 0, 1, for Ui 7^ 0, 1 and i chosen arbitrarily. 
Thin continuous lines correspond to power laws -J xja. 



NO (with A,, 7^ 1) become pairwise degenerated since the 
systems is split in two identical part by the rt^ = 1 plateau 
[Figs.|3;b),(c)]. Even in this case we find that the OPDM 
decays as a power law ~ ^^|a:;/a|~^/^ for large x in 
the regions away from n;, rij 7^ 0, 1 (Fig. 

The universal behavior of the OPDM at large-x above 
shows that although the ID HCB gas does not exhibit 
BEC in the TL quasi-long-range order is present 

and a large occupation of the lowest NO can occur. In 
the periodic case the NO are plane waves, so that their 
occupation and the MDF coincide. The results for the 
large-x behavior of the OPDM in the periodic case imply 
that in the TL n^j^g scales as \J Nh at constant p, and 
diverges as for fc ^ and A^j, 00, i.e., in the 

same way as in the case without a lattice 0, . 

In the trapped system, due to loss of translational in- 
variance, the NO and the MDF do not coincide. To 
obtain the behavior of the lowest NO in the TL, we 
study how it scales when the strength of the trap (or the 
number of particles) is changed keeping the characteris- 
tic density constant. In Figs. |2Ia)-(c) we show the re- 
sults obtained for the lowest NO in harmonic traps where 
the curvature of the confining potential was changed by 
one order of magnitude. It can be seen that a scaled 
NO (/j" ~ R^/^(tP exists, which docs not change when 
any parameter of the system is changed keeping p con- 
stant. Even when regions with Ui ~ 1 are present [Figs. 
I2b),(c)], where the two lowest NO are degenerate, the 
scaled NO exist. Here, the length scale C set by Va de- 
termines the length L = B'^C, over which a nonvanishing 
density is present. B'^ depends on the characteristic den- 
sity p and the power a of the confining potential. 
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FIG. 2: Scaled lowest NO and density profiles (thin dotted lines) for harmonic traps with (a) p = 0.55, Ai'i,=101 (continuous 
line), Afi,=30 (dashed line); (b),(c) p = 3.0, A^i,=551 (continuous line), A''j, = 167 (dashed line); (d) shows 10°! vs A'';,, in traps 
with power a = 8 of the confining potential, for p = 1.0 (-I-), 2.0 (x), 2.25 (V, a, rii — 1 region is present); the continuous lines 
correspond to power laws \^Nt. 



The scaling factor R is defined for any confining po- 
tential as i? = y/NhOja. A measure of the scaling of the 
NO can be obtained by studying the area below them. If 
the above-mentioned scaling is valid, since this area is ex- 
pected to depend only on /5, we have / dx\'^^ {x/C) |/C ^ 
R'^Ha/O E^\^^\=Sp implying that = ^^^A^ 

{Sp and S'p depend only on p). Figure |2Id) shows the re- 
sults obtained for J2i \ vs Ni, (at constant^ in traps 
with power a = 8 of the confining potential [l6j. It can 
be seen that already for iVb > 100 the expected power 
law \/ Nb is present, confirming our previous statements 
for large enough systems. 

With the findings above, the leading behavior of the 
lowest NO Ao = X^ij 4>^Pi34>^j can be evaluated in the TL 
and for a given p as follows. Replacing the sums in Ag 
by integrals (L 3> a) one obtains 



Ao ^ (l/a^) 



dx 



^ 0O(a:)A^0O(2/) 

dy-. 



iC/af'R-' 



\{x-y)/a\-y^ 



dX 



dY 



|A-y|-i/2 

(6) 



where we did the change of variables x = X(, y = YC,, 
and = The integral over X, Y depends only 

on the characteristic density. Then for a given confining 
potential with power a, and D'~ depend only on p, 
demonstrating that Aq scales in the TL as y^T^. The 
same analysis can be done with the MDF, where instead 



of normalizing it by the system size (as usual for homo- 
geneous systems) we normalize it by the length scale set 



-ife(i-j) 



Consid- 



by the potential [uk = {0'/C)J2ij ^ 
ering the large- a; form of the OPDM, and repeating the 
reasoning above, one obtains that nk=o also scales as 
or v^iVb, for constant p, in the TL. 

So far we have analyzed the scaling of Ao and 
«/c=o/\/C7^ in tlic TL. We discuss in what follows its 
relevance for finite size systems. In Fig. 13 we plot 
A' = X/ ^JC,/a for the first two NO (a) and n'^.^^ ~ 
n-k^o I a/C/o (b) as a function of the characteristic den- 
sity for two traps with v ery different confining poten- 
tials with power a = 8 [16j. As Fig. shows, finite 
size corrections to the leading behavior are very small 
since the values of A and nfc=o almost scale like in the 
TL. Actually, the results for the NO are indistinguish- 
able after the region with = 1 appears in the sys- 
tem, which is the point where the degeneracy sets in Fig. 
I2a)- Wc find the finite size corrections to A' and n'j,^g 
to be 1 / \fCja so that Aq / \f^Ja ^ Cp — Ep/ \f^Ja and 
rife=o/ yC7« ^Fp- G~pl^JZla, where Cp, Ep, Fp, and Gp 
depend only on p. In addition Ep and Gp can be positive 
or negative depending on the value of p. As an example 
we plot in the inset of Fig.|2Ib) 6Xq ~ Cp — \q/ \/C,/a for 
p = 1.0 {Ep is positive) and (5ri'^^Q = Fp — nk=o/ \/C/a for 
p ~ 2.25 (Gp is negative). Similar results were obtained 
for harmonic traps and the homogeneous system, in the 
latter case changing (/a by N. 
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FIG. 3: A' of the two lowest NO (a) and n'^^o (b) vs the characteristic density (p) for traps with power of the confining potential 
a = 8, Vga* = 2.0 x W'^^t (continuous line) and Vsa* = 1.0 x 10~^^t (dashed line). The inset in (b) shows SX'o for p = 1.0 (x) 
and (Jnjj^o for P = 2.25 (+) vs {C^/a)~^^^ (see text); the continuous lines show linear behavior. 



Finally, wc study the largc-77 behavior of the NO oc- 
cupations (A^). In contrast to the large- a; behavior of 
the density matrix, we do not find a universal feature 
in the large-77 behavior of for arbitrary values of the 
characteristic density p. However, for very small values 
of p we find that a universal power law develops in the 
large-77 region of A^, as is shown in Fig. 01 The power- 
law decay is in this case of the form A^ = AN^/ri"^, where 
An,, depends only on the number of particles in the trap 
independently of the confining potential, as is shown in 
Fig. ^ Since this occurs only for very low values of p 
we expect this behavior to be generic for the continuous 
limit. This, to our knowledge, has not been discussed 
before. In the latter limit the high momentum tail of the 
MDF was found to decay as Uk ~ 1^1"''^ for HCB in a 
harmonic trap 0|, and for the Lieb-Liniger gas of free 
and harmonically trapped bosons for all values of the in- 
teraction strength [l8|. Our results for the MDF (inset 
in Fig. 2)) show that the large-A: behavior of Uk for low p 
is also universal, irrespective of the confining potential. 

At this point it is important to remark that the univer- 
sal behavior and scaling relations shown in Figs. ^^01 ap- 
pear already at moderate number of particles, and hence, 
are relevant for experiments. However, the singu- 
larity in Tifc, well known from the homogeneous system, 
is difficult to see explicitly in such situations, in contrast 
to previous claims [l9j. Fitting power laws ~ for 
finite systems could lead to wrong conclusions about the 
large-x behavior of the OPDM, as shown in Fig. [3 For 
very low fillings (A^;,=ll), a "power-law" behavior with 
f3=l may be seen before the /?=4 is established for large 
k. Increasing the number of particles leads to a decrease 
of /3, up to (3=0.6 (A^b=401). Hence, the power /? depends 
strongly on the number of particles and cannot be under- 
stood as reflecting any universal property of the system. 
Power-law behavior disappears only when rn reaches 1 in 
parts of the system {Ni,=591 in Fig. [SJ. 
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FIG. 4: Occupation of the NO vs rj in systems with = 
10000: Nb = 21 periodic case (dashed-dotted line), TV;, = 11, 
V2a^ = 5.0 x IQ-^^t (dashed line), and Nt = 11, Vga** = 
5.0 X 10~'^*t (dotted line). The inset shows Uk vs ka for the 
same parameters and the same notation. Thin continuous 
lines correspond to power laws and in the inset. 



In summary, we have shown that quasi-long-range or- 
der is present in ID HCB on the presence of a lattice, 
with a universal power-law decay of the OPDM, inde- 
pendent of the power of the confining potential. Further- 
more, we have shown how the occupation of the lowest 
NO and the value of the MDF at zero momentum are 
determined by the large distance behavior of the OPDM. 
Even in the cases where a region with ni=l builds up 
in the middle of the system we find that both quantities 
scale proportionally to ^/T^ (at constant p). A further 
universal power-law decay has been found for the eigen- 
values of the OPDM (A,,) for large values of 77 at low 
densities (A,, ~ V~^)- It translates into a corresponding 
power-law decay of the MDF (n^ ~ 1^1"'*) at large mo- 
menta also independently of the power of the confining 
potential, pointing to scale invariancc in the ultraviolet 
limit of the continuous case. 
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FIG. 5: MDF for systems with V2a^ = 3.0 x 10"^t and Nb = 
11 (continuous line), A^i, = 101 (dashed line), and Nt = 401 
(dotted line), A^;, = 591 (dashed-dotted line). Accompanying 
thin continuous lines correspond to power laws: for Nt = 
11, k-" '' for Nt = 101, and for A^^ = 401. 

Note added. A HCB gas has been realized very recently 
on ID lattices by Paredes et al. [20| . 
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